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ABSTRACT 

We present high spatial resolution near-infrared spectra and images of the nucleus of Centaurus A 
(NGC 5128) obtained with NAOS-CONICA at the VLT. The adaptive optics corrected data have 
a spatial resolution of 0'.'06 (FWHM) in K- and O'.'ll in H-band, four times higher than previous 
studies. The mean velocities and velocity dispersions of the ionized gas ([Fell]) are mapped along 
four slit positions. The [Fell] emission line width decreases from a central value of a ~ 320 km s _1 
to a — 120 km s _1 at r ~ 0'.'5, and exceeds the mean rotation within this radial range. The observed 
gas motions suggest a kinematically hot disk which is orbiting a central object and is oriented nearly 
perpendicular to the nuclear jet. We model the central rotation and velocity dispersion curves of the 
[Fell] gas orbiting in the combined potential of the stellar mass and the (dominant) black hole. Our 
physically most plausible model, a dynamically hot and geometrically thin gas disk, yields a black 
hole mass of Mbh = 6.lto'g x 10 7 M Q . As the physical state of the gas is not well understood, we also 
consider two limiting cases: first a cold disk model, which completely neglects the velocity dispersion, 
but is in line with many earlier gas disk models; it yields an Mbh estimate that is almost two times 
lower. The other extreme case is to model a spherical gas distribution in hydrostatic equilibrium 
through Jeans equation. Compared to the hot disk model the best-fit black hole mass increases by a 
factor of 1.5. This wide mass range spanned by the limiting cases shows how important the gas physics 
is even for high resolution data. Our overall best-fitting black hole mass is a factor of 2-4 lower than 
previous measurements. A substantially lower Mbh estimate when using higher resolution kinematics 
was also found for many other black hole mass measurements as HST data became available. With 
our revised Mbh estimate, Cen A's offset from the Mbh-c relation is significantly reduced; it falls above 
this relation by a factor of ~ 2, which is close to the intrinsic scatter of this relation. 

Subject headings: black holes: general — galaxies: individual(NGC 5128) 



1. INTRODUCTION 

Galaxy merging, the formation of stellar spheroids, nu- 
clear star-formation and the fueling of nuclear black holes 
appear to be all linked, forming a central theme in build- 
ing galaxies. Most galaxies where these processes are 
currently acting, are so far away that the 'sphere of in- 
fluence' of the black hole cannot be spatially resolved and 
that detailed studies of stellar populations are impossi- 
ble. Yet, our own Galactic Center alone tells us that 
galaxy centers tend to become increasingly more inter- 
esting w hen observed at higher and higher spatia l resolu- 
tion fe.g. lSchodel et alJl200l iGenzel et al1l20f)! . When 
zooming into the nucleus of Centaurus A (NGC 5128), 
recent HST imaging (NICMOS) and ISAAC observa- 
tions have also revealed smaller and smaller 'sub -systems ' 
fScnreier et al ] IT^ lMa7cV)m et al .1120001 l2O0l . Cen A, 
the closest massive elliptical galaxy, the nearest recent 
merger, and one of the nearest galaxies with a signif- 
icantly active nucleus provides a unique laboratory to 
probe the interconnection between these phenomena on 
scales that are contained in the central resolution clc- 
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ment in any other object of its kind. 
The intricate dust lane that hides the center of Cen A 
has not allowed optical high-resolution spectroscopy with 
HST. Infrared (PR) spectroscopy is thus increasingly im- 
portant to open up the regime of dust shrouded nuclei 
and get accurate black hole mass measurements for these 
objects. 

Cen A is esp ecially interesti n g sinc e the black hole mass 
deduced by iMarconi et all l|200lD (M bh = 2.0±f° 4 x 
10 8 M(n) lies a factor of ten above the Mbh-c r elation 
ijFerrarese fc Merrittjl2000l iGebhardt et alJl2000|) . This 
is the largest offset from the relation measured to date. 
It is important to check this value, since it can help us 
to understand the coevolution of black holes and their 
surrounding bulges in more detail. 

We have initiated a program to stu dy the central p arsec 
of Cen A using N AOS/CONICA iRousset et alJ Il998t 
iLenzen et aTI IT998) at the Very Large Telescope (VLT). 
The adaptive-optics assisted imager and spectrograph 
provides us with near infrared data from 1 to 5 /xm at the 
diffraction limit of a 8m class telescope. The resolution 
is thus nearly fourfold that of HST in K-band. 
The distance to C en A is still u nder discussion. In a com- 
prehe nsive review llsraelllpQ^ gives a value of 3.40±0.15 
Mpc. iTonrv et all l)2001j) finds a value of 4.2±0.3 Mpc 
from I-band su rface brightness fluctuations, and recently, 
iReikubal 1)20041) derived a distance of 3.84±0.35 Mpc from 
Mira period-luminosity relation and the luminosity of the 
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tip of the red giant branch. Here, we assume a distance 
of D=3.5 Mpc to be consistent w i th the blac k hole mass 
measu rements of iMarconi et al.l l)2001|) and iSilee et aP 
( 2005) that also used that value. 

The paper is structured as follows: in Section 2 we 
present the observational strategy and the data reduc- 
tion. Section 3 describes the treatment of the spectral 
data. Section 4 presents the dynamical modeling and 
the results for the black hole mass for Centaurus A, and 
Section 5 discusses the implications of these results. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Adaptive Optics Observations 

Near infrared observations were performed in 2004 
March 28 and 31 with NAOS-CONICA (NACO) at the 
Yepun unit (UT4) of the very large telescope (VLT). 
NACO consists of the high-resoluti on near-infrared im- 
ager and spectrograph CONICA ijLenzen et alJ 11998(1 
and the Nasmyth Adaptive Optics System (NAOS) 
l)Rousset et al.lfl998|) . It provides adaptive-optics cor- 
rected observations in the range of 1-5 /jm with 14" to 
54" fields of view and 13 to 54 mas pixel scales. 
The data were taken in visitor mode and seeing during 
observations was in the range 0"3-0"8 (as measured by 
the seeing monitor in V-band), with clear/photometric 
conditions. 

Not seen in the visible, the active nucleus at the center 
of Centaurus A is a n unresolved source in K-band of 10.9 
mag as detected bv IMarconi et all (|2000). There are no 
potential reference stars bright enough (jtik <14 mag) 
for the wavefront correction at a distance of < 30" to the 
nucleus, necessary for a good quality of correction at the 
nucleus. Therefore, we directly guided on the nucleus 
itself using the unique IR wavefront sensor implemented 
in NAOS. This strategy provides us the best possible 
wavefront correction in the vicinity of the active galactic 
nucleus (AGN). In fact we reach the diffraction limit of 
the VLT in K-band of FWHM 0"057 and are not far off 
in H-band with O'.'ll. During the observations the atmo- 
spheric conditions were stable and the performance of the 
IR wavefront sensor (WFS) was steadily very good. For 
observations in H-band we used the K-dichroic, i.e. all 
the nuclear K-band light for wavefront correction. While 
observing in K-band itself the only possibility to achieve 
a good performance of the WFS was to send 90% of the 
light to NAOS and only 10% to CONICA (i.e. use the 
N90C10 dichroic). This increases the exposure times by 
a factor of 10 and made it effectively impossible to go for 
spectroscopy in K-band. 

2.2. K-Band Imaging 

For the imaging we chose the strategy to jitter the field 
on five positions at a separation of 4" and take a sky im- 
age at a dark position in the dust lane at a distance of 
^170" South-East of the nucleus; this cycle was repeated 
4 times. The on-chip exposure time was 120 s, yielding 
a total exposure time on the nucleus of 40 min. The re- 
sulting K-band image is shown in Figure Q 
The atmospheric conditions were stable and the seeing 
was 0"5 at the start and 0'.'8 at the end of the observa- 
tions. 

2.3. H-band spectroscopy 




Fig. 1. — K-band image taken with NACO. The nucleus is clearly 
visible as the unresolved source (FWHM=0'.'06). North is up and 
East is to the left. 

We took H-band spectr a at 4 different positio n angles 
and chose three similar to IMarconi et al.l l)2001|) (P.A.= 
-44?5, 32?5, and 82?5) in order to complement their 
findings with better spatial resolution. The fourth 
slit position (P.A.=70?5) was chosen as to observe the 
nucleus and a foreground star simultaneously to monitor 
the point spread function (PSF) on a stellar reference. 
See Figure for the positioning of the slits. 
The observations were obtained with the 0"086 slit and 
a grating with a wavelength range from A = 1.37/mi to 
1.84/im. The pixel scale was 0"054 leading to a slight 
under-sampling but shorter exposure times. The disper- 
sion was 7.0 A pixel -1 , yielding a resolution of R=1500 
in H. At a given position angle, the observations went 
as follows: first, the loop for adaptive optics correction 
was closed. When a good performance of the WFS was 
reached, an acquisition image in H-band was taken. 
The slit was centered on the prominent nuclear peak 
with an accuracy of ±0.2pixel (±0"011). The actual 
observations consisted of two sequences of exposures at 
4 positions along the slit. This was done in order to 
perform sky subtraction and to avoid detector defects. 
The effective slit length is 14". The exposure time per 
frame was 300 s, leading to a total exposure time of 
40 min per slit position angle. 

Data reduction was performed using standard IRAF 
routines. The frames were first bias corrected and flat- 
fielded with spectroscopic lamp- flats. Then cosmic rays 
were rejected and the frames were wavelength calibrated 
and corrected for distortions using spectroscopic arc 
lamps. 

To get the final 2-D spectrum the frames were aligned. 
This is done by shifting all frames to the same nuclear 
position. Finally, the frames were sorted by quality of 
wavefront correction, using both the width of the peak 
and the level of continuum flux. 



2.4. PSF -reconstruction 

A difficult but crucial part of adaptive optics obser- 
vations is the assessment of the point spread function 
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Fig. 2. — Slit positions overlaid onto the H-band acquisition 
image. North is up and East is to the left. From North to South 
the slit positions are 32?5, 70?5, 82?5, and -44? 5. The sign of 
the position along the slit used in the rotation curves is indicated 
by "+s" and "-s" . Note the foreground star that is in the slit at 
P. A. =70. 5° in the lower right. 

(PSF). The PSF is highly dependent on the quality of 
the wavefront correction, quantified e.g. by the Strehl 
ratio. The Strehl ratio is set by the observable proper- 
ties of the reference object (flux, size, contrast) but also 
by the atmospheric conditions (seeing, coherence time); 
it is therefore changing with time and needs to be mon- 
itored throughout the observations. In case of the K- 
band image a separate PSF reference star was observed 
directly after the nucleus with the same WFS setup. This 
star was chosen fr om the 2MASS point source catalogue 
ijCutri et alJl2003|) to match Cen A's nucleus as closely 
as possible: in angular proximity, magnitude and color. 

Since the acquisition for NACO observations takes a 
non-negligible amount of time, going back and forth to 
the PSF star is a tedious and time-consuming task. On 
the other hand the measurement of a separate PSF-star 
at a given difference in time and position on the sky does 
not guarantee to give a good approximation on the actual 
on-source PSF. We therefore came to the conclusion that 
it is most suitable to measure the PSF in the science 
frame on the unresolved nucleus. 

For emission line spectroscopy the strategy is to have 
the nucleus in the slit at all position angles and measure 
the PSF on this unresolved point-source for each frame. 
We tested this approach by choosing one position angle 
(P.A.= 70°) such that the nucleus plus a foreground star 
are simultaneously observed, and indeed, the widths of 
both light profiles are similar (compare Table [T]l . 

We describe the normalized PSF empirically by a sum 
of two Gaussians; one narrow component describing the 
corrected/almost diffraction limited PSF core (<tdl) and 
one broader component which we later attribute to the 
seeing halo (er s ): 

PSF (r) = -Ve- 2 /^ L + (l^e-TV**, 

DL 27r(7 s 

where F is the ratio of the flux of the narrow component 
and the total flux of the PSF (F = fluxDL/flux to tai)- The 
quantity F provides a rough approximation of the Strehl 



ratio (S) which gives the quality of an optical system; S is 
defined as the observed peak flux divided by the theoret- 
ically expected peak flux of the Airy disk for the optical 
system (S = peakfluxDL/peakfluxAiry)- For the following 
analysis it is sufficient to measure the quantity F, which 
also gives an estimate of the quality of the adaptive op- 
tics correction. 

The width of the fitted broad component can be com- 
pared with the seeing that is given in the header for 
each frame, as measured by the seeing monitor in V- 
band. One has to adjust the seeing estimates to the 
same wavelength, as the resolution depends on wave- 
length A, as 6{X) ~ A -1 / 5 . The individual components 
are given in Tableland are compared to the header in- 
formation (adapted to H-band) . The agreement between 
<7 S and ctdim is satisfactory in all cases. Notice also the 
good agreement between the width of the nucleus and 
the star at P.A.=70?5. The adaptive optics correction is 
optimised for the position of the nucleus and is worse at 
the position of the star due to anisoplanetism. 
Moreover it is obvious that the width of the narrow com- 
ponent depends on the seeing conditions. 
Figure 0] shows the integrated flux over the two- 
component model PSF shown in Figure Note that 
50% of the flux lie within a radius of O'.'l and 90% of the 
flux within 0"3. 

In the PSF model we did not account for the undersam- 
pling of the observed PSF, since the sampling problem is 
negligible compared to the general uncertainty in adap- 
tive optics observations. 



TABLE 1 
Reconstruction of the PSF 



Object 


PA. 


FWHMdl 


FWHM S 


FWHM DIM(H) 


F 


Nuc 


32? 5 


O'.'ll 


0"37 


0'.'37 


0.23 


Nuc 


-44? 5 


0'.'15 


tf'48 


0'.'48 


0.18 


Nuc 


82? 5 


0'.'15 


0'.'53 


0'.'49 


0.15 


Nuc 


70? 5 


O'.'ll 


0^'34 


0'.'35 


0.22 


Star 


70? 5 


0'.'12 


0"36 


0'.'35 


0.31 



3. RESULTS 

3.1. Nuclear spectrum 

A nuclear H-band spectrum centered on the continuum 
peak (Fig. EJ and extracted from a 0'.'054 x 0'.'086 aper- 
ture is presented in Figure [5J The spectrum exhibits a 
power-law continuum with three [Fell] lines: Ail.534/im, 
A2l.644/im, and A3l.677/^m. We use the strongest line 
(A2l.644/mi) for our kinematical studies. 

Linc-of-sight velocities, FWHMs, and surface bright- 
nesses along each slit were obtained by fitting single 
Gaussians to the [Fell] emission line in each row of the 
continuum-subtracted two-dimensional spectra. The cor- 
responding slit positions are shown in Figure[5]and listed 
in Table [T] 

3.2. Gas Kinematics 

We measure the gas kinematics on the ionised [Fell] 
line AI.644/im out to around ±0'.'6. Single Gaussians 
provide a good fit to the emission lines and are used to 
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Fig. 3. — Decomposition of the PSF for the slit position angle 
82? 5. The data are over plotted by the 2-component PSF-model 
(solid line), where the dashed line indicates the corrected peak and 
the dotted line represents the uncorrected seeing halo. The bottom 
panel shows the residual between model and data. 




r [arcsec] 



Fig. 4. — Integrated flux over the model PSF is shown for the 
slit position angle 82? 5 (to be compared with Figure 151 . 




1.0 .... i .... i .... i .... i ... , 
1.50 1.55 1.60 1.65 1.70 1.75 

Fig. 5. — Nuclear spectrum in H-band extracted from an aperture 
of 0^054 X 0'.'086 at a slit position of 32? 5. The three [Fell] emis- 
sion lines are indicated at the rest-frame wavelengths: Ail.534/im, 
A2l.644/im, and A3l.677/im. 
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Fig. 6.— Observed [Fell] A 2 1.644/mi line on the NACO H- 
band spectrum at a slit position of -44? 5. To make the gas velocity 
curve more visible, each column of the spectrum was divided by the 
maximum value of the emission line in that column. The spectrum 
was then resampled on a 3 times finer grid by means of bilinear 
interpolation. The line indicates the systemic velocity of 532 km 
s — , which is shifted to zero for this plot. 

measure the position and width of the line. The fit is per- 
formed in IDL using a non-linear least squares fit to the 
line and the errors are the 1-a error estimates of the fit 
parameters. The center of the continuum peak is taken 
as a reference for the systemic velocity. We find a value 
of v sys =532 km s ; in good agreement with the value 
v sys =532±5km s _1 measured by iMarconi et alJ (2001) 
from their gas kinematical data. 

For the central 0'.'3 only the four highest quality frames 
are considered in order to make use of the full resolution. 
The quality of the frames, i.e. their correction quality, is 
estimated on the basis of the width and the peak value 
of their continuum peak. Outside 0'.'3, when the lincflux 
drops, all frames are taken into account and three pixels 
are binned to enhance the signal-to-noise ratio. 
The rotation of the gas can be directly seen in Figure 
where the line intensity was normalised by the peak in- 
tensity in each column. The rotation curves, velocity dis- 
persion curves and emission-line surface brightness pro- 
files are shown in Figured for all slit positions. 

The velocity dispersion is directly measured as the 
width of the lines. We corrected for instrumental broad- 
ening of 65 km s , measured from the skylines. 
The excellent spatial resolution of the NACO data 
is demonstrated in Figure |S1 where we compare the 
NACO data points t o the kinematical data published by 
IMarconi et all(|2001|) and lSike et al.l(|2005|) . The NACO 
velocity dispersions are often considerably smaller than 
the ISAAC velocity dispersions measured at the same lo- 
cation. The information given in lMarconi et al.l l)2001|) is 
not sufficient to find the reason for this discrepancy and 
we believe that our data points are correct. 

3.3. The Emission-Line Surface Brightness 

The intrinsic surface brightness distribution of emis- 
sion lines £ is an important ingredient in the model 
computations because it is the weight in the averaging 
of the observed quantities. It would be ideal to have an 
emission line image with higher spatial resolution than 
the spectra. Unfortunately , the resolution of th e HST 
[Fell] narrow band image ( Mar coni et alJ l2000|) is not 
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Fig. 7. — Rotation curves (left panel), velocity dispersion profiles (middle panel) and emission line surface brightness profiles (right panel) 
along the 4 slit positions. For the position angles -44? 5 and 70? 5 some of the points shown in the surface brightness profile are missing for 
the v and <r profiles. 

convolve them with the PSF of the observations. The 
two-dimensional gas distribution is assumed to be a cir- 
cular disk. We find that the intrinsic surface brightness 
of a disk is well fitted by a double exponential profile 

- r/ro J- Iie- r/ri 
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Fig. 8. — Comparison between the spatial resolution of the 
NAC O data (filled diam onds), the ISAAC data (open trian- 
gles) (Marconi ct al. 2001), and the GNIRS data (open squares) 
ISilge et aH2003) . The NACO and ISAAC data correspond to gas 
kinematics while the GNIRS data to stellar kinematics. In the left 
panel the rotational velocity is shown and in the right the velocity 
dispersion. 

good enough to mimic the NACO data and therefore 
our approach is to match the emission-line fluxes as ob- 
served in the spectra. We extracted the emission line 
surface brightness from [Fell] directly from the spectra 
along the four slit positions (see Figure [7| right panels). 
In order to get a parametrisation for the intrinsic sur- 
face brightness, we test different functional forms and 



where ro and rj are the scale radii, and Iq and Ii are the 
scale factors for the t wo components. This parametri- 
sation was also used in lMarconi et alJ (|2001^l . Since the 
observed emission line surface brightness along the slit 
at different position angles depends on the inclination of 
the circular gas disk (with respect to the observer) as well 
as on the PSF we fit the intrinsic emission-line fluxes for 
a median disk inclination of 45° given the observational 
setup and conditions. Th e projected major ax is of the 
[Fell] disk is at P.A.=32?5 ijMarconi et alJ200C(l and gets 
the highest weight in the fit, since the surface brightness 
along the major axis stays unchanged when the disk is in- 
clined. The parameters that after PSF convolution best 
fit the surface brightness distribution along the differ- 
ent slit positions are Ii/I =0.018, r =0'/02, and ri=0"22. 
This parametrisation for the disk surface brightness is 
used in the dynamical models at all inclinat i on a ngles. 
Compared to the fit values in iMarconi et al.1 <|200lD our 
scale radii are smaller due to higher spatial resolution. 
Given the shape o f the [Fell] emission from the HST 
narrowband image ijMarconi et al J 12000(1 . the axis ratio 
of a possible gas disk is ~1:2 which corresponds to an 
inclination angle of 60°, given a circular disk structure. 
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However, as found before bv lMarconi et all 1)200 If) there 
is a mismatch between the photometric major axis of 
the [Fell] gas disk and its kinematical major axis (line-of- 
nodes), which they find to lie between —10° > ( > —14°. 

4. DYNAMICAL MODEL 

In this section we outline the models describing the mo- 
tions of the ionised gas in the combined gravitational po- 
tential of the (putative) central black hole and the stars. 
This modeling also requires a careful discussion of the 
gas geometry, as well as the physical state of the ionised 
gas, whose measured velocity dispersion far exceeds the 
expected thermal broadening. 

4.1. The Stellar Mass Model 

The first step in the dynamical modeling is to estimate 
the stellar contribution to the central potential from the 
stellar surface density. For the deprojection of the ob- 
served surface brightness distribution into the stellar lu- 
minosity density we applied the Multi-Gaussian expan- 
sion (MGE) method lIMonnet. Bacon, fc E msellcm fll992t 
lEmsellem. Monnet. &: Baconlll994|) . Specifically, we ob- 
tained an MGE fit to a composite set of 2-dim ensional K- 
band images using the method and software o flCappellaril 
( 2002). The MGE fit was performed using the NACO K- 
band image (0" < r < 7"), the NICMOS F222M image 
ISchreier et alJll99l (2" < r < 16"), an d the 2MASS 
Large galaxy atlas (LGA) K-band image ijJarrett et al.l 
2003) for larger radii (8" < r < 365"). 
The sky-subtraction was performed relative to the 
2MASS LGA image, which is taken to be the sky-clean 
reference. The flux calibration of the NACO and NIC- 
MOS data is also referenced to the 2MASS image. 

To get a parametrisation of the NACO K-band 
PSF we fit a full MGE model to the composite im- 
age of NGC 5128 (without giving the code a prior 
PSF estimate). We then selected the central four 
Gaussian components, that clearly described the un- 
resolved nucleus, as the PSF. Its parametrisation is 
PSF(R) = ^!liGiexp[-R 2 /(2cr i * 2 )]/(27rcr i * 2 ), and the 
numerical values of the relative weights Gi (normalised 

such that J2i=±pi = 1) an d of the dispersions a* are 
given in Table El This PSF was then fixed in the MGE 
fit to the composite image, where we neglected PSF 
convolution of the NICMOS and 2MASS images since 
these are only used at radii r > 2" and the NACO image 
has a much higher resolution. 



TABLE 2 

MGE PSF PARAMETRISATION 

i Gj a* FWHM* 

1 0.0004 0'.'005 0'.'012 

2 0.1590 07031 0'.'074 

3 0.3034 07062 07145 

4 0.5372 0'.'143 0'.'336 



Figure shows a comparison between the observed 
photometry and the MGE model along four different 
position angles in the galaxy, while Table gives the 
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Fig. 9.— Comparison between the NACO (0" < R < 7"), 
NICMOS (2" < R. < 16"), and 2MASS LGA (8" < R < 365") 
K-band surface brightness (open squares) of Centaurus A and the 
corresponding best-fit MGE model (solid line) as a function of ra- 
dius R. On the left the MGE fit is shown for inclination angles 
between 0° and 90° , while the right panels show the corresponding 
residuals. 

corresponding numerical values of the analytically de- 
convolved MGE parametrisation of the galaxy surface 
brightness 

E(x',y') = V L \ exp --^ (x' 2 + ^] , 

where (x , y ) are the coordinates on the plane of the 
sky and N is the number of the adopted Gaussian com- 
ponents, having total luminosity Li , dispersion <7j, and 
observed axial ratio 0.8< q i <1.0. In Table the central 
unresolved component (which we attribute to the AGN) 
is removed. We only account for the stellar light distri- 
bution in the Multi-Gaussian expansion model. 

TABLE 3 

MGE PARAMETRISATION OF THE DECONVOLVED K-BAND SURFACE 
BRIGHTNESS OF CEN A 



i 


U [xlO 9 L , K ] 


o~i 


% 


1 


0.113 


0"45 


0.882 


2 


0.297 


T/48 


0.812 


3 


1.13 


3727 


1.000 


4 


1.25 


6" 68 


1.000 


5 


5.07 


14"3 


0.879 


6 


15.5 


3870 


1.000 


7 


19.6 


7876 


1.000 


8 


23.7 


13079 


0.807 


q 


. 5.24 


365"8 


0.800 



The total luminosity hi, the dispersion a* , and the axial ratio q ; 
for the nine Gaussian components that are needed to fit the surface 
brightness of Cen A in the Multi-Gaussian expansion model. We 
overlap three K-band images (NACO (0" < R < 7"), NICMOS 
(2" < R < 16"), and 2MASS LGA (8" < R < 365")) to model 
the surface brightness from 0" to 365" (see Fig. |5J. The central 
component seen in Fig.|5]is omitted here. 

Using this radial profile of K-band (volume) emissivity, 
we fix the stellar mass-to-light ratio by constructing an 
isotropic spherical Jeans model and matching it over 
the radial range 3" < r < 10" to the stellar kinematics 
(\A>* +o-l) published recently by iSike et all (|2005|) . 
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The innermost 3" are excluded from the fit to minimise 
the influence of the black hole. We find a best-fitting cen- 
tral mass-to-light ratio M/L K = (0.72 ± 0.0 4) M^/L^ 
in agr eement with the best-fitting model of ISilee et alJ 
( 2005) . This agreement of the mass-to- light ratio from 
Schwarzschild and Jeans mode l ing is supported by the 
recent work of ICappellari et all l|2006j) . 



4.2. Geometry and Kinematics of the [Fell] gas 

The central [Fell] velocity curves (Figure |7J) suggest 
that we see gas rotating in a flattened geometry around 
a central mass concentration, presumably dominated by 
a black hole. The velocity dispersion in the very center 
is quite high (~300 km s _1 in the central pixel), well 
in excess of the mean rotation seen. In part this high 
dispersion may be attributed to rapid but spatially un- 
resolved rotation. 

As the underlying physical origin of the observed high 
gas dispersion is not clear a priori, we consider at least 
three geometric models: 

1. the gas lies in a geometrically thin, kinematically 
cold disk and is in Keplerian motion around the 
black hole, 

2. the gas forms a geometrically thin, kinematically 
hot disk that is in radial hydrostatic equilibrium, 

3. the gas lies in a spherical distribution of collision- 
less cloudlets and can be modeled through Jeans 
equation. 

Note that cases 1. and 3. are extreme physical assump- 
tions and we consider them as limiting cases. In all three 
cases the gas is moving in the combined potential of 
the surrounding stars and of the central black hole; the 
self gravity of the gas is negligible (given the estimated 
mass of the ionised gas is ~ 1O 3 M l|Marconi et alJ 20011 
and the mass of the innermost Gaussian representing the 
stars is already ~8x 1Q 7 M Q (see Table EJ). The under- 
lying stellar mass distribution is the same in all models, 
with M/L fixed by the stellar kinematics at larger radii 
(see previous section). 

The mismatch between the kinematical a nd photometric 
major axis of the gas disk, found before by Marconi et al. 
(2001) introduces a degeneracy between the inclination 
angle (i) of the gas disk and the position of its line-of- 
nodes (C). We fear that our long-slit data will not be 
sufficient to constrain both i and £. We therefore make 
use of the inclination angle of the radio jet which sets 
the lower limit of the disk inclination in the standard 
pictur e of an orthogonal disk-jet geometry. iTingav et alJ 
(1998) derived a value of 50° < i < 80°. Given this prior 
information, we set up our thin disk models. 

4.3. Model 1: Thin cold disk model 

We follow t he w idely used approach (e.g. 
iMacchetto et al.l ^997) which assumes the gas to 
lie in a thin disk around the black hole, moving on cir- 
cular orbits; its observed velocity dispersion is assumed 
to be solely due to rotation. We con structed a model 
using the IDL software developed in ICappellari et al.l 
(2002J. This modeling can deal with multiple component 
PSFs, different PSFs for different data-subsets, and 
a general gas surface brightness distribution. Pixel 
binning and slit effects are taken into account to 



generate a two-dimensional model spectrum with the 
same pixel scale as the observations. Like in the data 
analysis, the rotational velocity and velocity dispersion 
is determined by fitting simple Gaussians to each row of 
the model spectrum. We make no velocity offset correc- 
tion ( van der Marel et all Il997t iMacieiewski fc Binnevl 
l200lt IBarth et alJl2001^ . since the slit width (0"086) is 
comparable to the FWHM of the PSF. 
The predicted velocity and velocity dispersion profiles 
of this model depend on the intrinsic surface brightness 
distribution of the emission lines, the PSF, plus the 
following parameters: 

1. the inclination, i, of the gas disk (i = 0° is face-on, 
90° edge-on) 

2. the angle between the projected major axis of the 
disk (line of nodes), £, and the slit positions 

3. the black hole mass, Mbh; the stellar mass profile 
is fixed. 

We describe the "intrinsic velocity dispersion" of the gas 
in the model by a double exponential parametrisation of 
the form 



o-r = <? e 



-r/r a 



(1 



-r/ri * 



without using it in the dynamics. The parameters a®, e, 
7"o, and r\ are fixed through comparison of the observed 
dispersion profile to a flux-weighted convolution of the 
intrinsic velocity dispersion profile with the PSF and 
the size of the aperture. The parameters that give a 
good fit to the velocity dispersion at all slit positions 
in the cold disk model are, cto=140 km s" 1 , e ~ 0.25, 
ro = 0'.'9, and n = 0'.'2. We do not know the source for 
this high velocity dispersion and simply ignore it in the 
kinematically cold disk model. To reduce the number of 
free parameters in the model we evaluate the rotation 
curves for three fixed inclinations (45°, 60°, and 70°) 
and fixed the intrinsic surface brightness density of the 
gas disk beforehand (cf. Section l3~3^l . 
In Table 0] we give the best-fitting black hole masses 
for each inclination angle individually; we also give the 
resulting % 2 values to get the overall best-fitting model. 
The la, 2a, and 3cr contours are shown in Figure 1101 
for the two degrees of freedom, black hole mass and 
line-of-nodes. The three sets of contours correspond to 
the relevant inclination angles of the gas disk (i=45°,60°, 
and 70°) where the best-fitting inclination angle, i=45°, 
is indicated by a solid line. 



TABLE 4 
Best Fit Values - Cold Disk Model 



i M bh [M e ] 



45° (4.0 ±0.2) x 10 7 -27° ±4° 46.1 
60° (3.2 ±0.2) x 10 7 -28° ±4° 47.1 
70° (3.2 ±0.1) x 10 7 -36° ±2° 47.1 



The overall best-fitting cold gas model is obtained at 
an inclination angle of 45°. Figure ITTI shows the observa- 
tional data in comparison to the model. The parametri- 
sation of the surface brightness, shown in the third row, 
is fixed beforehand and not included in the fit here (com- 
pare Section 13. 3|) . The parameters that lead to this fit 
are C = -27° and M bh = 4.0 x 10 7 M Q . 
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-24 - 
-26 - 
-28 - 
-30 - 
-32 - 
-34 - 

-36 - 

-38 L 

2.0x10 7 




i = 70* 



2.5x10 7 3.0x107 3.5x10 7 

M B „ [M s ] 




= 45° 



4.0X10 7 4.5x10 7 5.0X10 7 



situation reads ijBinnev fc TremainelllQS'il Eq.4-64a): 



Fig. 10. — Fitting a cold gas disk: the x 2 contours are shown for the two parameter fit of black hole mass vs. projected major axis £ 
for the three inclination angles. The best cold model is at an inclination angle of 45° and the corresponding contours are plotted in solid 
lines. The best-fitting values are £ = -27° ± 4° and M bh = (4.0 ± 0.2) X 10 7 M Q . The three contours give the formal 68.3%, 95.4%, and 
99.73% (thick) confidence levels. 

However, the minimum \ 2 values are comparable for all 
three cases, i.e. their differences are smaller than the la 
level, and we therefore cannot constrain the inclination 
angle of the gas disk with our long-slit data. We conclude 
that the black hole mass for a dynamically cold, geomet- 
rically thin disk model is M b h = (^O^g) x 1O 7 M . 



R ®M) , R d* 



dR 



dR 



4.4. Model 2: Thin hot disk model 

This model is identical to the previous one except that 
we now interpret the "intrinsic velocity dispersion" as a 
gas pressure component. As mentioned above, we do not 
know the physical origin of this high velocity dispersion 
of the gas, nevertheless, we include it in the dynamical 
analysis. 

Any pressure support will require less rotation for dy- 
namical equilibrium at a given black hole mass. In other 
words, the black hole mass needs to be larger to cause the 
same mean rotational velocity. The classic approach to 
account for velocity dispersion is to ap ply an asymmetric 
drift correction Ce.g lBarth et al.1 1)2001[) '). but the approx- 
imate equations are only applicable if a/v is small. How- 
ever, in the case of Cen A, a/v even exceeds unity and 
we therefore chose another approach: we assume the gas 
disk to be geometrically flat but with an isotropic pres- 
sure; on this basis we construct an axisymmetric Jeans 
model in hydrostatic equilibrium to model the observa- 
tions. 

We assume that the [Fell] surface brightness of the gas 
disk T, g (x,y) reflects the tracer gas density p g (x,y), and 
that the gas moves in the joint potential <E> of the stars 
and the central black hole. The Jeans equation for this 



where R is the projected radius, and an and v$ are the 
radial velocity dispersion and the azimuthal velocity of 
the gas, respectively, and both are functions of R. 
We again parametrise the gas dispersion profile by a dou- 
ble exponential function of the form: 



~ r / r 0< 



(T R = a e ''"'(l + ee^ /ri ), 

and we fit for the best set of parameters (<7o = 140 km 
s-\ e ~ 0.25, r = 0'.'9, and r x = 0'.'2) to get the 
intrinsic dispersion profile. 

The assumption that the gas disk is infinitesimally 
thin does most probably not reflect the real physical 
properties but was chosen to eliminate line-of-sight 
integrations through a 3-dimensional gas distribution. 
Table [5] summarises the best fit parameters and the 
corresponding x 2 values for the different inclination 
angles of the "hot" gas disk model. The lcr, 2a, and 3a 
contours are plotted in Figure fHfl 

The best-fitting hot disk model favors a black hole 
mass of M bh = (6.1 ± 0.3) x 10 7 M© and £ = -27° ± 3°, 
at an inclination angle of 45° . It is shown in Figure IT21 
for comparison. The increased total \ 2 value (compared 
to the cold disk model) can be explained by the reduced 
number of degrees of freedom, since in the hot disk 
model the velocity dispersion and the rotational velocity 
are coupled through Jeans equation. However, the 
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Fig. 11. — The observational data is shown in comparison to the best fitting cold disk model at 45° inclination angle. The best-fitting 
parameters are £ = —27° and M^h = 4.0 X 10 7 Mq. The velocity dispersion as well as the surface brightness are parametrised by a double 
exponential function which were fixed beforehand and not fitted here. 



TABLE 5 
Best Fit Values - Hot Disk Model 



C 



i M bh [M Q ] 

45° (6.1 ±0.3) x 10 7 -27° ±3° 52.1 

60° (5.6 ±0.3) x 10 7 -25° ±2° 54.0 

70° (6.4 ±0.3) x 10 7 -25° ±2° 53.7 



difference between the two fits is small (< la confidence 
level) and the large velocity dispersion lead us to 
consider this model a better description of the physical 
properties at the center of Centaurus A. 



4.5. Model 3: Spherical Jeans model 

In this final modeling approach, we account for the 
fact that the high gas velocity dispersion is not consis- 
tent with the assumption of a thin disk. Alternatively, we 
assume the gas spherically symmetric distributed in in- 
dividual clouds that move ballistically (i.e. t co n>>td yn )- 
Given the observational facts (rotation cu rves along cer- 
tain s lit positions, gas disk observed by iMarconi et al.1 
( 2000)), this model is not very likely to describe the phys- 
ical properties of the central region in Cen A. However, 
it gives an upper limit on the black hole mass, compared 
with disk models at non-extreme inclinations, and is im- 
plemented fairly easily. We construct a spherical Jeans 
model where we assume the following: 

1. the gas density is given by its emissivity 

2. the gas cloud distribution is spherical 



3. the gas clouds move in the potential given by the 
stars and the possible central black hole 

4. the stellar mass-to-light-ratio is constant through- 
out the relevant range (< 2") 

5. the stars are in spherical symmetry. 

Here, we construct a model with a Multi-Gaussian- 
expansion both for the stellar photometry and for the 
gas distribution. The gas surface b rightness is given by 
the HST [Fell] narrow band image (|Marconi et alJl200o|l 
and fitte d by an MGE model. 

Following Trcmainc et alJ l)1994l) in the spherical case the 
solution of Jeans equation for the projected rms velocity 
reduces to 

2G f°° p g (r)M tot (r) 



(v z +ct%(R) 



T£ ff (i?) 



-(r 2 -R 2 ) 1/2 dr, 



where T is the stellar mass-to-light ratio, S s (i?) is the 
surface brightness of the gas, p g (r) is the density of the 
gas and M tot (r) is the mass of the stellar body and the 
central dark object. 

Again, the stellar mass-to-light ratio is fixed at a value 
of T = O.72M /L , that we derived in section IP The 
best fitting black hole mass that we find with the spher- 
ical Jeans model is (1.0 ± 0.5) x 10 8 M Q . The compari- 
son between the measured and the modelled rms velocity 
(V v 2 + cr 2 ) is shown in Figure ITTl 

5. DISCUSSION 

The availability of adaptive optics instrumentation 
opens up new realms of ground based observations. 
Drawing on kinematic data with unprecedented spatial 
resolution we present a dynamical model of the black 
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Fig. 12. — The best-fitting hot disk model with parameters £ = —27° and M^h = 6.1 X 10 7 Mq. In this model, v and a are interwined 
through Jeans equation. Note that despite the higher black hole mass (compared to the cold disk model), the model curves drop beyond 
0'.'3 and match the observed data better. The model of the surface brightness is fixed beforehand and not fitted here. 
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Fig. 13. — Fitting the 'hot' gas disk: the \ 2 contours are shown for the two parameter fit of black hole mass vs. projected major axis C, 
for the three inclination angles. The best hot model is at an inclination angle of 45° and the corresponding contours are plotted in solid 
lines. The, best fitting values Jire C = -27°.± 3^,and M b K =f (6.1 ± 0.3) X, ipjM©..The three contours give the formal 68.3%, 95..4%_and 

g^TO^tififi) EcM8fficcOfcvki™ aurus A. lire spatial res- olution that we reacn m our H-band spectra is 0.11 
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Fig. 14. — Comparison between the measured rms velocity (filled 
diamonds) and the spherical Jeans model (open diamonds) for the 
gas. The best fitting black hole mass is M bh = (1.0±0.5) X 10 s M . 

(= 1.8 pc at a distance of 3.5 Mpc) that is a factor 
of 3 to 4 higher than from t he previous ground b ased 
observations with ISAAC bv iMarconi et all {200 At 
large radii (> 0'/3) our data are in agreement with 
their gas kinematical data, but our high resolution data 
reveals the steep gradient of the rotation curve in much 
more detail (see e.g. Figure 0. The rotation curves 
show a smooth behavior and the slit position that shows 
the fastest rotation (P.A.=— 44?5) is nearly perpendic- 
ular to the jet direction (P.A.~ 51°). This suggests 
the picture of gas orbiting in a flattened geometry 
(presumably a disk) around the central black hole, with 
the jet pointed along the disk angular momentum vector. 



5.1. Black hole mass 

In this paper we present three different dynamical 
modeling approaches to describe the gas kinematics 
in the central - 2" ~ 33pc of NGC5128. Two of the 
models are clearly conceptually inconsistent with the 
data, either neglecting the dominant velocity dispersion 
of the gas (Model 1; cold disk), or the disk-like geometry 
of the velocity field (Model 3; isotropic, spherical Jeans). 
We have laid out these extreme models as limiting 
cases, which would result in a broad mass range, 
2.7 x 10 7 M© < M b h < 1.5 x 10 8 M Q , and to demonstrate 
how important the gas physics is, even in light of high 
resolution data. Our physically most plausible model 
is a hot disk model at an inclination angle of 45°; the 
corresponding black hole mass is Mbh = 6.1+jj'g x 10 7 Af Q 
with a line-of-nodes at £ = —27° ± 4°. We consider 
this model the best-fitting since it accounts for the high 
velocity dispersion of the gas combined with a disk-like 
gas structure, as suggested by the smooth rotation 
curves. The assumption of a geometrically thin disk is 
not physically motivated but was chosen to eliminate 
line-of-sight integrations through an underconstrained 
3-dimensional gas distribution. 

Our disk model meshes well with other constraints on 
the ce ntral geometry of NGC 5128: e.g. iTingav et alJ 
(1998) derived a jet position angle of ~ 51°. If we as- 
sume that the accretion disk and surrounding gas disk is 



at right angles, we would expect its line-of-nodes to be 
at ~ —39°. Moreover, from the direction of the jet and 
counter-jet they derived a value for the jet inclination 
of 50° < i < 80° with respect to the line-of-sight. Our 
best-fitting inclination angle of 45° is close to their lower 
value, but our long-slit data do not tightly constrain the 
inclination angle. This was also the main source of un- 
certainty for the blac k hole mass in the previous study by 
IMarconi et al.l 1(20011) . If we were to consider a gas disk 
inclination to be as face-on as 25° (inconsistent with the 
jet inclination), the black hole mass would increase to 
- 1.2 x 10 8 Af© in the hot disk case. 

The unsettled question of inclination will presumably 
be solved with the analysis of integral field spectroscopy 
data taken with SINFONI at the VLT. The full 2-D ve- 
locity field data provide an excellent means of modelling 
black hole masses from gas as well as stellar kinematics. 

Our best-fit value for the line-of-nodes is only ~ 12° 
away from the "expected" val ue of ~ —39° , muc h 
closer than the value derived bv IMarconi et all |2001), 
-10° > C > -14°. Nonetheless, this confirms that the 
kinematical line-of-nodes does not coincide with the 
projected major axis of the gas di sk (P.A.~ 33°) seen 
in Paa and [Felll w ith NICMOS l|Schreier et alJll99Sl 
IMarconi et al.ll2000|) . 

The value that we derived for the black hole mass 
from the best-fitting hot disk model (6.llg;| x 10 7 M Q ) 
is significantly lower than the values bot h from 
stellar kinematics presented by ISilge et all (|2005f) 
(M bh = l-8±oj x 10 8 M Q for i =45°) and the previ- 
ous gas kinematical study of IMarconi et alJ 1(200 ID 
(M b h = 2-0±l l x 10 8 M Q ) obtained at 3-4 times lower 
resolution. This decrease in black hole mass estimate 
with higher resolution data is in line with the decrease 
in derived black hole masses when HST data became 
available. Our limiting case model of a spherical gas dis- 
tribution modeled through Jeans equation gives a black 
hole mass of Mbh = ( 1-0 ± 0-5) x 10 8 M(p which agree s 
both with ISilge etafl {2005) and IMarconi et all fe OOl). 
Howev er, in their gas dynamical modeF lMlircoiii et alJ 
l(2001() assumed a cold thin disk and did not account for 
the pressure support of the gas; assuming a cold disk 
model, we find the best-fitting black hole mass to be 
M b h = (4.0i°;„) x 10 7 M©, which is almost a factor of 7 
lower than their value. 



5.2. Relation of black hole mass versus galaxy properties 

This confirmation of a fairly high black hole mass 
compared to a fairly low stellar velo city dispersion 
of cr* = 138 km s _1 ( Silg e et aP |2005() quantitatively 
reduces but qualitatively confirms the offset of Cen- 
taurus A from the Mhh-ff re lation (Ferrarese fc Merrittl 
120001: IGebhardt et~alT hoOOj) . The black hole mass 
predicted by this relation would be around 3 x 10 7 Af© 
and our best-fitting black hole mass lies a factor of 
~ 2 above this. This offset is close to the observed 
scatter of the Mhh-g which is a factor of 1.5 (0 . 2dex) 
((Ferrarese fc Merrittl 1200(1 IGebhardt et all 12000). In 
the case of the lowest black hole mass supported by our 
cold disk model (M bh = (4.0±?:§) x 10 7 M Q ) the black 
hole falls nicely onto this relation. 
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To compare the black hole mass to the bulge mass of 
Centaurus A, we applied a spherical Jeans model to 
the whole galaxy as seen in K-band, u s ing th e stellar 
kinematics (v» and ct») of iSike et all l)2005j) out to 
~ 100" and derived a total mass of the spheroid of 
M sp h = (2.5 ± 1) x 10 11 M Q . This mass is in agreement 
with the value M sph = (4 ± 1) x 10 n M w de rived by 
iHui et all l|1993|) and iMathieu fc Deionghd (|1999|) . 
Given this spheroid mass, Cen A lies b elow the M^- 
Meuigc relation (e.g. Marine k, Rixl 12004) (i.e. it has a 
fairly low black hole mass compared to its bulge mass) 
but is not a striking outlier in this relation. Taken 
together, it seems that Cen A foremost has a very 
high M sp h/<7» ratio among ellipticals (which implies a 
very low concentration), rather than being an outlier in 
the relations to Mbh- This low concentration may be 
explaine d by the fac t that Cen A is known to be a z~0 
merger (|Israellll998l) . 

This paper demonstrates that near-IR adaptive optics 
instrumentation provide excellent data and make it pos- 



sible to explore the central regions of dust enshrouded 
galaxies. With a spatial resolution of ~ 0'.'06 in K-band 
and ~ O'.'ll in H-band we are now able to resolve the 
radius of influence of black holes even in more distant 
galaxies from the ground. 
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